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INTRODUCTION 


The  problem  of  determining  the  state  of  stability  of  a slender, 
flexible  missile,  modeled  mathematically  as  a free-free  elastic  beam 
subjected  to  a thrust,  has  been  the  subject  of  several  recent  inves- 
tigations [1]  - [9],  It  is  of  particular  interest  that  Wu  [9]  has 
determined  that,  in  the  absence  of  a feedback  directional  control, 
the  mechanism  of  instability  of  the  missile  system  described  in 
references  [2],  [4],  and  [8]  is  not  flutter  but  rather  divergence. 

This  discrepancy  appears  to  result  from  the  form  of  approximate 
solution  assumed  in  those  investigations.  For  example,  to  solve 
approximately  the  fundamental  nori-sel f-adjoint  boundary  value  problem 

y'v  + Q(xy')1  - w2y  = 0,  0<x<l , (1) 

y"  = y"‘  =0  at  x = 0,1,  (2) 

where  y1  = dy/dx,  etc.,  and  Q and  w denote  the  dimensionless  thrust 
and  frequency  parameters,  respectively,  Beal  [2]  employed  the 
Galerkin  procedure.  He  assumed  a solution  that  included  the  rigid 
body  motions  of  translation  and  rotation,  and,  as  a result,  always 
obtained  two  zero  frequencies.  However,  it  is  very  easy  to  verify 
that  equations  (1)  and  (2)  admit  a rigid  body  mode  of  rotation  when- 
ever Q>0,  i.e.,  there  exists  only  a single  zero  frequency.  Subsequently, 
Wu  [9]  reported  that  the  loss  of  stability  of  the  nature  of  a 
structural  mode  that  degenerates  to  the  rigid  body  mode  of  rotation 
when  Q ► 0 rather  than  by  the  onset  of  flutter  which  ensues  upon  the 
coalescence  of  two  non-null  natural  frequencies  at  a positive  critical 
value  of  the  thrust  parameter  Q. 


The  objective  of  the  present,  investigation  is  to  examine  a 
related  missile  problem  from  a somewhat  simpler  viewpoint,  so  as  to 
expose  rather  clearly  through  elementary  mathematical  relationships 
the  qualitative  stability  characteristics  of  a flexible  missile.  The 
system  will  be  modeled  mathematically  as  a system  comprising  three 
rigid,  weightless  bars  joined  by  elastic  springs,  carrying  a con- 
centrated mass  at  the  center  of  each  bar,  and  subjected  to  a cir- 
culatory thrust,  whose  direction  may  be  adjusted  through  a variation 
in  the  value  of  a so-called  tangency  coefficient. 

2.  THE  EQUATIONS  OF  MOTION 

Consider  the  system  depicted  in  Figure  1 that  consists  of  three 
rigid,  weightless  bars  of  equal  length  Si  with  elastically  hinged 
joints  that  exert  the  following  linear  restoring  moments:  c(<J>^  - i}^) 
and  c(4>  - <f>  ),  where  the  <f>  - ' s , j = 1,2,3,  are  the  ( assumed  small) 

J c J 

angles  formed  between  the  horizontal  and  each  of  the  respective  bars 

and  c denotes  a constant  stiffness  parameter.  The  neutral  equilibrium 

configuration  corresponds  to  <p.  = 0.  The  concentrated  masses  m.  are 

J J 

affixed  at  the  center  of  each  bar,  and  x and  y are  the  coordinates 
of  the  position  of  the  central  mass  m^  relative  to  the  fixed  xy  - 
coordinate  frame.  A thrust  of  constant  magnitude  Tq  with  its  line  of 
action  specified  by  the  tangency  coefficient  a,  whose  vector  form  F 
is  F = Tq  (cos  + sin  a^j),  where  i and  j are  unit  vectors  parallel 
to  the  fixed  x - and  y - axes,  respectively,  is  applied  at  the  left 
free  end  of  the  system.  It  is  assumed  that  the  system  can  translate 
in  the  vertical  y - as  well  as  the  horizontal  x - directions. 
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In  the  usual  manner,  one  can  derive  the  following  expressions  for 
the  kinetic  energy  T,  the  potential  energy  V,  and  the  generalized 


forces  Qx,  . . . , Q3: 


T “ 2"  (ml  + m2  + m3^x2  + 2 ^ml  + m2  + m3^2  + 8fT'l5'2^l  + 

+ 1 (m1  + m 3)?,2^  + 1 ny2<5>2  + \ V*  si"  + 

+ 2 ^1  ~ m3^X  Sln  <f>2<V  2 "V*  Sin  ^3  " 

- ny  y cos  <f>^  - j (m^  - iy  cos  + 

+ 1 nyy  cos  4>34>3  + ~ny2  cos  - 4^  )<y2  + 

+ ~ mi2  cos  (<J>2  - 4)3^2^>3  ’ ^ 

V = ~ C - 2<t>1  <f>2  + 24>|  + 4)3  - 2<ty3)  , (4) 

Qx  = T0  cos  “*i»  = Tq  sin  aV 

= T0H  sin  (1  - a)^  , Q2  = ]rV  Sin  (<J>2  - ct^),  Q3  * 0.(5) 
In  view  of  the  forms  of  equations  (3)  - (5),  the  Lagrange  equations 


may  be  written  as 


S $ ‘ Q*  • Tt  (f  ■ V 

.1  fil  ) . -3i-  + 9V_  = Q.  , j = 1 ,2,3 

dt  3<j)j  J 
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Therefore,  substituting  equations  (3)  - (5)  into  equations  (6),  one 
finds  that  the  linearized  equations  of  motion  for  the  given  system  are 

(m-|  + m2  + m3)x  = T0  (7) 

(m1  + m2  + m3)y  ~ (m-,  - m3)£$2  + ~ m3&$3  = «T0<t»1  , (8) 

4 + ~ rn^Ax^  - j m^y  + ~ m^2^  + c (^  - <*>2)  = T0£(l  - a)^  , 

1 (m]  + m3)^2^2  + 2 " rV£*<t2  " ^mi  ' m3^  + 4 V^l  + 

1 1 (9) 

+ - m3A2$3  - c<{)1  + 2c<j>2  - cd>3  = — Tq  ^(d>2  - <*♦] ) » 

4"  m3£2^3  ' 2 m3£X(J,3  + \ m3£y  + \ m3e,2^2  + c^3  ' C<J>2  = 0 » 

where  the  familiar  small  angle  approximations  have  been  made  and  higher 
order  terms  have  been  neglected. 

Solving  equations  (7)  and  (8)  for  x and  y,  one  finds 

x = pTo  (10) 

and 

y = u[~  m1A?1  + “ (rn1  - m3)£$2  - “ m3£$3  + aTQ  ^ ] , (11) 

where  y = l/(m^  + m2  + m3).  Insertion  of  equations  (10)  and  (11)  into 
equation  (9)  yields  the  following  set  of  three  equations  of  motion: 

-■  pn»i  (m2  + m3)Jl2$1  + [c  - -yTo£(l  - a)(m1  + 2m?  + ?m3)]4>1  + 

+ — ym^  (m2  + 2m3)£2$2  - c<+>2  + ym^m3£2<f3  = 0 , 

i-ym^£2(m2  + ?m3^1  “ Cc  “ jr  Va(m2  + ^3)^  + 
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+ 1 pJL2 (m-jm2  + 4m'-jm3  + m2m3)c|>2  + [2c  - ^ 

- - yTo£(m2  + 2m3 ) ]4>2  + ~ ym3£,2(2m-|  + m2)$3  - = 0 , 

4 Vm-)m3S'2$-]  + ym3aTQ  ~ ym3(2m^  + m^i2'^  ~ d<} >2  + 

+ 1 ym3(m1  + m2)£2$3  + (c  - ^ ym^A)^  = 0 . 

It  is  convenient  to  use  a dimensionless  form  of  equations  (12). 
For  this  purpose,  the  following  definitions  may  be  made: 

= -j  m 9 ni^  — m, 

t = ot,  o2  ~ m£2/4c(l  + y-j  + y3)  » 

Q = q/(l  + y1  + y3)  , q = TQ  &/2c  . 

It  is  now  a straightforward  process  to  verify  that  the  system  of 
equations  (12)  assumes  the  form 
3 

l [Amn^n  + (Cmn  + QDmn)4>n]  = 0,  m = 1,2,3 
n=l 

where  now  $n  = d24)n/d-r2,  and 
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= y i ( 1 

f * 

A 12 

= P1  (1 

+ 2y3) , 

A13  ~ 

21 

= y i ( 1 

+ 2y3) , 

A22 

= yl  + 

4^3  + y3 

A = 
, h23 

‘31 

= *Y3 

t 

A32 

= ^3° 

+ 2y ^ ) , 

A33" 

11 

= 1, 

C12 

= - 1, 

C1 3 = 

21 

= - 1, 

C22 

= 2, 

C23  = 

31 

= 0, 

C32 

= - 1, 

C33  = 

Dn 

= (a  - 1 ) (2  + 

P1  + 2y3), 

D12 

-0,  d13  = 0, 

°21 

= ot  ( 1 + ?P3). 

°22 

= -(l+2y3),  D23  = 0, 

°31 

= a^3 

°32 

= 0,  D33  = - y3  . 

(15) 

THE 

EREQUENCY  EQUATION 

If 

a solution  of 

the  form 

<Dn(t)  = Xne 

io)T 

i 

n = 1 ,2,3 

(16) 

where  the  Xn's  are  constants,  is  substituted  into  equation  (14),  then 
a system  of  homogeneous  algebraic  equations  in  the  X n ' s is  obtained. 
This  system  has  a non-trivial  solution  if  and  only  if  the  determinant 
of  the  coefficient  matrix  vanishes.  Expansion  of  this  determinant 
yields 

Pou>4  - p2w2  + P4  = 0,  (17) 

where 

p = 8y  y (1  + p + y ) - y.y  Q[3  + 4(y  + 2y  ) + (y  + y ) (y  + 5y  ) - 

o 13  13  13  I 3 131  3 

- a(l  + y1  +■  P3)2]. 

P2  = 4 (y i + 4y1y;J  + y3)  - Q[3y ^ + 1 ly 3 + 36y1y3  + y]  + 13y2  + 

+ 34y^y2  + 10y|y3  - a(t ^ + H M 3 + + 22'J1,J3  + + 

(IS) 

+ 10y2y3  + lOy^y2)]  + 2y3Q2  [1  + 3y^  + 4y ^ f 1 1 p 3 + 

+ y2  + 3y 2 + 71j11j3  + 3y2y3  - <*0  + 21^  + 4^3  + 

+ 7yjU3  + yj  + 3y2  + 3y^'3  + 3y|y3)], 
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P4  = - Q(1  - a)[p3(2  + v1  + 6y g + -jU 3 + 4y2)Q2  - (2  + + 

+ 1 1 V! 3 + 4p -j VI3  + 10y2)Q  + 3 + + 5)^]  . 

A priori,  one  would  have  expected  the  frequency  equation  (17)  to 

be  a bi -cubic  polynomial  in  m2.  However,  one  can  easily  verify  that 
the  inertia  matrix  A (Amn)  is  singular  which,  therefore,  accounts 
for  the  fact  that  the  frequency  equation  is  a biquadratic.  Consequently, 
the  system  under  consideration  is  a quasi  - dynamic  system  [10],  which 
is  characterized  by  the  existence  of  an  internal  constraint.  To 
determine  this  constraint,  one  merely  adds  the  equations  for  m = 1 
and  3 obtained  from  equation  (14)  and  then  subtracts  the  equation  for 
m = 2 from  the  result  to  obtain 

[2  - Q < ( 1 + y ) - Q ( 1 - a)(2  + u]  + 2y3)]<J>1  - [4  - Q(1  + + 

+ (2  - P3Q)^3  = 0 . (19) 

It  is  obvious  from  equation  (18)  that  p^  = 0 when  Q - 0.  Thus, 

the  frequency  equation  (14)  leads  to  w2  = 0,  p*/p*  , where 

<-  o 

p*  = 8y1y3(l  + y]  + P3),  P*  = 4(u1  + 4y1y3  + y3)  . 

The  value  w2  = 0 corresponds  to  a rigid  body  rotation  of  the  system. 

4 • STABILITY  CONSIDERATIONS 

If  Q has  an  arbitrarily  small  positive  value,  then  p^  f 0 and 
0 is  no  longer  a root  of  the  frequency  equation  (17).  Thus,  one 
■ist  roni'lude  that  the  system  can  no  longer  undergo  a rigid  body 
10‘ition.  As  the  value  of  Q is  increased  above  Q = 0,  an  eigencurve 
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emanating  from  the  origin  in  the  Qw2-plane  must  enter  either  the  first 
or  second  quadrant,  which  implies  that  the  system  is  either  stable  or 
unstable,  respectively.  Clearly,  if  one  can  establish  a relationship 
between  Q and  w2  that  is  invalid  in  a neighborhood  of  the  origin  in 
the  Qw2-plane,  then  some  conclusions  regarding  the  stability  of  the 
system  can  be  drawn. 

To  accomplish  this  for  small  thrusts  Q,  the  frequency  parameter 
w2  is  expanded  as  a power  series  in  the  perturbation  parameter  Q as 
fol 1 ows : 

co2  = Qa-j  + Q2a2  + ( Q 3 ) , (20) 


where  the  coefficients  , c^,...  are  to  be  determined.  Substituting 
equation  (20)  into  the  frequency  equation  (17)  and  following  the 
steps  of  the  classical  perturbation  process,  one  finds 


Q(1  - a) (3  + + 5y3) 

4(p1  + 4Mip3  + p3) 


provided  that  0 < Q <<  1.  From  this  equation  it  is  evident  that  the 
sign  of  w2  is  determined  solely  by  the  magnitude  of  the  tangency 
coefficient  a.  Specifically,  if  a < 1,  then  w2  < 0 and  the  system  is 
divergent  for  an  arbitrarily  small  positive  thrust.  Hence,  the  value 
of  the  critical  thrust  of  divergence  is  = 0 when  a < 1.  If  a > 1 
(a  super-tangential  thrust),  then  w2  > 0 and  the  system  is  stable  for 
sufficiently  small,  positive  thrusts. 

If  the  thrust  is  a tangential  force  (a  = 1),  then,  provided  that 
Q > 0,  p^  again  vanishes,  as  is  evident  from  equation  (18).  In  this 
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f 


(22) 


case,  the  frequency  equation  yields 

u)2  = 0,  p2/po  , 

where 

p0  = 2u1p3[4(1  + + y3)  - Q(1  + y i + 3y3  + 2y-|U3  + 2u|)]  , 

P2  = 4(y1  + 4y -jy 3 + y3)  - 2Q[y-|(l  + 7y3)  + y|(l  + 12y-|)]  + 

+ 2yiy30  + 2y3)zQ2.  (23) 

A change  in  the  sign  of  w2,  as  given  in  equation  (22),  will  be  signalled 

by  the  vanishing  of  either  p2  or  pQ.  Specifically,  if  P2  = then, 

according  to  equation  (23)  one  has 

y^O  + 2m3)2Q^  - [^0  + 7y3)  + P20  + 12y])Qb  + 

+ 2(y^  + 4y.jy3  + Vj)  = 0 » 

whereas,  if  p =0  the  critical  thrust  for  divergence  is  given  by 

o 

4(1  + y-j  + U3) 

°b  = 1 + u,  + 3k3  * V3  * ?l,3  ' <2E> 

It  must  be  noted  that  p^  can  vanish  whenever 

y3(2  + + 6y3  + 2y^y3  + 4y2)Q2  - (2  + ^ + llp3  + ^3  + 

+ 10y2)Qb  + 3 + y^  + 5y3  = 0 , (26) 

which  also  provides  a critical  thrust  for  divergence.  The  values  of 

Q obtained  from  equation  (26)  are  completely  independent  of  the 
b 

tangency  coefficient  a. 
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The  condition  for  the  determination  of  the  critical  thrust  for 


I 


I 


i 


flutter,  namely  that  of  coalescence  of  a pair  of  roots  of  the 
frequency  equation  (17),  is  p|  - 4^  = 0 which,  by  virtue  of 

equation  (18),  can  be  expressed  as  the  following  quartic  polynomial 
in  Q: 


5 

l 

n=l 


r C"-1 
n e 


= 0, 


(27) 


where  the  r 's  are  lengthy  polynomials  in  p ,y  , and  a which  will  not 

I J 

be  recorded  here. 


5.  SOME  EIGENCURVES 

To  identify  regions  of  stability,  flutter,  and  divergence  in  the 
stability  maps  that  will  be  presented  later,  it  is  first  efficacious  to 
solve  the  frequency  equation  (17)  and  then  to  examine  some  typical 
eigencurves,  such  as  those  shown  in  Figures  2 to  4.  These  figures  were 
prepared  for  three  values  of  the  tangency  coefficient:  a = 3/10,  4/5, 

and  3/2,  in  the  case  of  p-|  = = 1 • In  agreement  with  the  conclusion 

drawn  from  equation  (21),  one  sees  from  Figures  2 and  3 that  the  branch 
of  the  eigencurve  emanating  from  the  point  (co2,Q)  = (0,0)  extends 
immediately  into  the  second  quadrant.  Thus,  the  system  is  unstable  by 
divergence  for  arbitrarily  small,  positive  values  of  Q.  Indeed,  in  the 
case  of  a = 3/10,  the  value  of  the  critical  thrust  of  divergence  must 
be  Qb  = 0.  In  the  case  of  a = 4/5  (Figure  3),  one  has  m2  = 0 at  Q = 0, 
a)2  < 0 for  0 < Q < 0.4127,  and  oj2  > 0 for  0.4127  < Q < 0.506.  The 
bounds  in  these  inequalities  are  obtained  from  equation  (26). 
Consequently,  one  may  conclude  that  the  system  diverges  whenever  the 
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thrust  parameter  Q is  in  the  interval  0 < Q < 0.4127  and  that  it  is 
stable  when  Q is  in  the  interval  0.4127  < Q < 0.506.  At  Q = 0.506, 
the  frequencies  of  the  first  and  second  modes  merge,  which  signals  the 
onset  of  instability  by  flutter.  For  Q > 1.2,  flutter  gives  way  to 
di vergence. 

In  Figure  4 plotted  for  a = 3/2,  the  branch  of  the  eigencurve 
initiating  at  the  origin  enters  the  first  quadrant,  reaches  a 
maximum  value  of  w2,  and  then  decreases  until  it  intersects  the  load 
axis  at  = 0.4127. 

The  dashed  horizontal  lines  appearing  in  Figures  2 to  4 are 

horizontal  asymptotes  which  represent  the  value  of  Q at  which  w2 

becomes  infinite.  This  value  of  Q,  hereafter  designated  as  Qa,  is 

determined  from  the  condition  p =0,  which,  by  virtue  of  equation  (18), 

o 

leads  to 

8(1  + y-j  + y3) 

q — . 

a 3 + 4(y1  + 2y3)  + (y-j  + y3)(y-j  + 5y3)  - a(l  + y-j  + y3)2 

6.  STABILITY  MAPS 

Having  in  mind  the  observations  made  in  the  preceding  section,  one 
now  can  identify  the  regions  of  stability,  divergence,  and  flutter  in 
qa-plane?  In  Figures  5 to  7,  stability  maps  are  shown  in  the  qa-plane 

*The  thrust  parameter  q was  defined  in  equation  (13).  It  is  more 
desirable  to  consider  the  qa-plane  rather  than  the  Qa-plane  since  the 
quantity  Q tends  to  mask  the  influence  of  the  mass  distribution  of  the 
system  on  the  value  of  the  critical  thrust. 
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over  the  range  0 * a < 2 for  (a),  = 5,  y = 1,  (b)  y-j  = y^  = 1 , and 

(c)  y^  = 1 , y = 5,  respectively.  The  regions  of  stability,  flutter, 

and  divergence  are  labelled  with  the  symbols  SR,  FR,  and  DR. 

It  should  be  noted  th?t . for  a < 1,  the  minimum  critical  thrust  of 

divergence  is  q,  = 0,  whereas,  for  a > 1,  its  value  is  positive  and  is 

b 

associated  with  the  smallest  of  the  two  roots  of  equation  (26).  For 
anti-tangential  (a  < 0)  and  super-tangential  (a  > 1)  thrusts,  flutter 
does  not  occur. 

From  Figures  5 to  7,  it  is  clear  that  there  exists  a value  of  a, 
say  a*,  at  which  qb  = q&.  The  value  of  a*  may  be  computed  from  the 
requirement  that  equation  (27)  and  equation  (26),  written  now  in  the 
form 

bl  + t>2  Qb  + b3  Qj  = ° » 

where 

bi  = 3 + y-|  + 5e3#  b2  = - (2  + y-j  + 1 1 y 3 + 4y^^  + lOy^), 

b3  = y3(2  + y-|  + 6y3  + 2yj  -j  y3  + 4y3), 

possess  a common  root  in  the  thrust  parameter.  This  requirement  leads, 
according  to  Sylvester's  dialytic  method  of  elimination  [11],  to  the 
condi tion 
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r5 

0 

b3 

0 

0 

0 


0 

b. 


0 

0 


0 

0 

b. 


0 


= 0 


Expansion  of  this  determinant  produces  a polynomial  in  a from  v/hich  the 
value  of  a*  may  be  determined  numerically  for  given  values  of  p ^ and 

V 

Consequently,  one  may  now  observe  that  for  a sub-tangential 

thrust  (a  * 1)  for  which  0 < a < rx*.  flutter  cannot  occur  and  the 

minimum  critical  thrust  of  divergence  is  q^  = 0.  However,  for  a 

sub-tangential  thrust  for  which  a*  < a < 1 > flutter  can  occur  (e.g., 

see  Tigure  3).  For  a in  this  latter  range,  the  system  is  first 

unstable  by  divergence  for  0 < q < , where  q^  = Qbl  (1  + + Vj), 

with  Q denoting  the  smallest  root  of  equation  (26),  and  it.  becomes 
bl 

stable  whenever  q.  < q < q . where  q^  = Qo  (1  + P,  + p ),  with 

u-j  e-j  e-|  » j 

Q denoting  the  smallest  real  root  of  equation  (27).  For  q < q 
e1  y el 

< q , the  system  is  prone  to  flutter,  and  for  q < q it  is  divergent 
e2  2 
once  again. 

The  stability  maps  in  Figures  5 to  7 indicate  that  a stable 
flight  at  relatively  low  values  of  the  thrust  parameter  q can  be  achieved 
provided  that  the  direction  of  the  thrust  is  controlled  in  such  a way 
that  a > 1.  The  critical  super-tangential  thrust  is,  of  course,  qb  . 
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Suppose  next  that  the  value  of  a is  only  slightly  greater  than 

unity.  If,  as  the  value  of  q is  increased  to  , the  value  of  a is 

changed  to  a value  very  slightly  less  than  unity,  then  a stable 

motion  with  thrusts  greater  than  q would  be  possible.  In  this 

bl 

second  region  of  stability,  the  maximum  permissible  sub-tangential 

* * 

thrust  would  be  q < qg  < q^  , where  the  values  of  q0  and  q^  can 

be  obtained  from  equations  (27)  and  (24),  respectively.  For  at  only 

★ _ 
slightly  less  than  unity,  qg  ~ q^  . For  example,  with  = y^  = ' > 

one  finds  qL  = 1.2381  and  qf  = 3.  Thus,  a significant  increase  in 
bl  bl 

the  admissible  magnitude  of  the  thrust  appears  possible,  provided, 
certainly,  that  the  control  system  can  accomplish  the  necessary  reduction 
in  the  value  of  a in  a period  of  time  that  is  sufficiently  brief  so  as 
to  avoid  large  angular  changes  relative  to  the  equilibrium  configuration 
of  the  system. 

Finally,  it  is  of  interest  to  consider  the  influence  of  the 

distribution  of  mass  of  the  system  on  the  value  of  the  critical  super- 

tangential  thrust.  To  do  this,  one  may  first  presume  that  the  total 

mass  of  the  system  is  assigned  some  value.  Then,  in  terms  of  the 

dimensionless  mass  parameters,  one  may  write  1 + = 8 + 1 , 

where  3 is  assigned  a fixed  value.  Therefore,  y^  - 8 - , where 

0 < y3  < B.  Substitution  of  y^  = 3 - into  equation  (26)  and 

introduction  of  Q.  = q./(B  + 1)  lead  to 
b b 


-1 
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P3[2  + B + (5  + 2B)p3  + 2vi|lq§  - (B  + 1 )[2  + B + 2(5  + 2b)p3 

+ 6y*]qb  + (B  + l)2(3  + B + 4^)  = 0 . (28) 

In  Figure  8,  the  variation  of  qb  , i.e.,  the  lowest  root  of 
equation  (28),  versus  y3  for  B = 1,2,..., 5 has  been  plotted.  All  the 
curves  show  essentially  the  same  behavior,  namely,  the  value  of  q^ 
decreases  monotonical ly  as  p3  increases,  with  the  difference 
A = q (g1)  - qb( b" ) , where  B'  < 3"  and  y3  held  fixed,  increasing  as 
the  difference  B"  - B'  is  increased.  Consequently,  significantly 
greater  critical  thrusts  result  when  the  tail  portion  of  the  system 
is  more  massive  than  the  head  portion.  The  lowest  critical  thrusts 
are  obtained  when  the  head  is  more  massive  than  the  tail. 
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